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The paper deals with the ground and the first excited state of the polaron in the one dimensional 
Holstein model. Various variational methods are used to investigate both the weak coupling and 
strong coupling case, as well as the crossover regime between them. Two of the methods, which are 
presented here for the first time, introduce interesting elements to the understanding of the nature 
of the polaron. Reliable numerical evidence is found that, in the strong coupling regime, the ground 
and the first excited state of the self-trapped polaron are well described within the adiabatic limit. 
The lattice vibration modes associated with the self-trapped polarons are analyzed in detail, and 
the frequency softening of the vibration mode at the central site of the small polaron is estimated. 
It is shown that the first excited state of the system in the strong coupling regime corresponds to 
the excitation of the soft phonon mode within the polaron. In the crossover regime, the ground and 
the first excited state of the system can be approximated by the anticrossing of the self-trapped and 
the delocalized polaron state. In this way, the connection between the behavior of the ground and 
the first excited state is qualitatively explained. 



I. INTRODUCTION 

Ever since LandauB suggested that the electron can 
be trapped by the deformable lattice, strongly coupled 
electron-phonon systems have been the subject of inten- 
sive examination. Besides the investigations of those sys- 
tems in which the lattice is coupled to the whole elec- 
tronic band, there has been significant interest in the 
physics of a single polaron, in which the electron and 
the associated lattice deformation form a quasi-particle, 
spatially and spectrally decoupled from the rest of the 
system. Lattice degrees of freedom make even a single 
polaron problem a many-body one. The analytical and 
numerical examinations of most electron-phonon models 
are thus difficult. For this reason, even the simple Hol- 
stein modela (suggested in 1959) is still being investigated 
in recent works. Various methods have been proposed in 
order to calculate the polaron ground state of the Hol- 
stein model. Almost exact results (except for the adia- 
batic limit) have been obtained with the quantum Monte 
Carlo calculations,™ the global-local method,0 the den- 
sity matrix renormalization g roup . method,Q and some 
exact diagonalization methods .r 2 lr'i 

The main goal of this paper is to determine the el- 
ements important for the qualitative description of the 
polarons in the whole range of electron-phonon coupling. 
Some of them, although already known, are found to be 
better understood when supplemented with additional 
details. 

The ground state of the Holstein system changes from 
the delocalized polaron state, in which the electron is 
nearly free, to the small and self-trapped polaron state, 
as the electron-phonon coupling g increases. These two 
opposite limits are usually identified as the weak and 
strong coupling regime, respectively. In the weak cou- 
pling regime the influence of the small lattice deforma- 
tion on polaron dynamics is very small, which makes the 
energies of polaron and electron hopping to neighboring 
sites similar, t po i < t. 

The exact ground state is an eigenstate of the system 



momentum, regardless of the couplingE Therefore the 
polaron ground state is delocalized for all parameters. 
However, in the strong coupling regime t po i becomes neg- 
ligible, leading to self-trapped polaron states. According 
to Ref. H the dynamics of the small self-trapped polaron 
can be separated into two time scales. On the short time 
scale, the lattice deformation is centered at some lattice 
site, and the electron can virtually hop among the neigh- 
boring lattice sites. Only after a certain number of such 
events (of the order t po i/t), the polaron as a whole tun- 
nels to a new central site. The localized polaron states 
thus may lead to a very accurate estimation of the po- 
laron ground state energy in the strong coupling regime. 
Accordingly, for the self-trapped polarons, an exact di- 
agonalization method of calculating the localized polaron 
states, rather than the translationally invariant ones, can 
be used. By comparison to the results of other methods, 
it is shown that this approach introduces only minor er- 
rors in the ground state energy of the self-trapped po- 
laron. Moreover, the localized polaron functions permit, 
unlike the translationally invariant ones, a separate anal- 
ysis of the electron and phonon properties of the polaron. 
The local electron density, the mean lattice deformation, 
and the on-site zero point motions of the self-trapped 
polarons can be calculated in this way. 

In the crossover (intermediate) regime, which is be- 
tween the weak and the strong coupling regime, no known 
perturbation calculation converges, which complicates 
the discussion of the polaron nature. Although it has 
been proved that [the change of the polaron ground state 
with g is smooth,EI the physics of the rapid crossover be- 
tween two opposite limits of the electron-phonon cou- 
pling, in the small interval of g's, is not completely clear. 
In Ref. |23| it has been claimed that the phonon excitation 
associated with the first excited state of the system is un- 
correlated to the electron in the weak coupling regime, 
while it is confined to the electron in the strong coupling 
regime. The transition occurs in the crossover regime 
in which, in addition, the energy difference between the 
ground and excited states becomes small. In the present 
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paper, the excited polaron states are treated by the new 
method which uses variational approach in order to de- 
fine and solve the generalized eigenvalue problem. Even 
if this method does not converge systematically, it does 
provide new results concerning the nature of the polaron 
ground and first excited state. 



II. GENERAL 



The Holstcin Hamiltonian reads 



H = ~t^24i,s (Cn+l,s + Cn-\,s) + TtoJ b^bn 



71,8 

9^2 c i,s c n,s (bl + b n ) 

n.s 



(1) 



It describes the tight-binding electrons in the nearest- 
neighbor approximation, coupled to one branch of dis- 
persionless optical phonons. s is the creation operator 
for the electron of spin s at lattice site n, and b\ is the 
creation operator for the phonon. t is the transfer (hop- 
ping) integral of the electron. Tiuj and g are the phonon 
and the electron-phonon coupling energies, respectively. 
The Holstein Hamiltonian depends only on two ratios 
of relevant energy parameters: g/fiio and t/%u), i.e., the 
results will use Tiuj as the energy unit. It is often con- 
venient to express the lattice vibrations in terms of the 
nuclei space and momentum coordinates, 

x n = x (bl + b n ) , p n = ip (bl - b n ) , (2) 



where xq = y/h/2Mu> and po = ^JMTiuj/2 are space 
and momentum uncertainties of the harmonic oscillator 
ground state. M is the mass of a nucleus and K denotes 
the spring constant, to 2 — k/M. By noting that the 
electron-lattice displacement coupling constant a, in g = 
axo, is independent of M, an alternative set of Holstein 
Hamiltonian parameters can be introduced, 



t — H 2 /2m e ia 2 



M 



a 
2k 



£_ 



which is convenient for the discussion of the adiabatic 
regime M ^> m e i . Here, m e i is the electron effective mass 
and a is the lattice constant. It is worth noting that t 
and E p (the binding energy of the polaron for t = 0) are 
independent of M and thus they are the only parameters 
relevant in the adiabatic limit. 

By using standard conventions for c\ s and b q , 

n n 

Eq. (Q) can be rewritten in momentum space and sepa- 
rated in two mutually commuting parts, H = Hk,q^o + 

Hq=0, 



Hk^o = -2t cos (ka) n kyS + hu ^ n q 

k,s q^O 

- g/VN cl +q , s c k Abq + bl q ) 

k,s,q^0 



H q=0 = huj h q=0 - gN e i/v N (b q=0 + b 1 ) 



where h kyS = <4, S C M' = b l b qi and N ei = Sfc.s^fe.s- 

The part which involves only the q = phonon mode, 
H q= o, can easily be transformed into the diagonal form, 
by using the unitary operator S q (£), 



S q (0 = exp , 
5-i (04=o4=o(0 = huh q=0 -(N el g) 2 /NhuJ, 



where £ = N e ig/vNhuj. Any eigenstate of H q=0 , 
S q= o(l;)\n q = ) , has the same mean total lattice deforma- 
tion, Xtot = J2n~Xn> 



xtot = ^oV^V(n g=0 |5' 9 io(0( & g =o + b q =o)S q= o(£)\n q=0 ) 
= x Q VN(C + = 2x N el g/huj , (3) 

which is independent of t. Since the q = part of the 
Hamiltonian commutes with H^^o, it can be concluded 
that Eq. (||) is also valid for all eigenstates of the total 
Holstein Hamiltonian. The phonon part of these eigen- 
states can be represented in the form of a direct prod- 
uct of two groups of states, the first one includes q ^ 
phonon modes, while the second includes only the q = 
phonon mode. This is useful because one can always 
check approximate computations by calculating Xtot, or 
include this property in the computation itself. This spe- 
cific property of the q — mode is not restricted to the 
particular dimension of the system, nor to the number of 
electrons. Moreover, it can also be found in some other 
models in which the electron-phonon coupling consists of 
the lattice deformation linearly coupled to the local elec- 
tron density. A hint in this direction was reported for 
the first time in Ref. 

The total momentum of the system, K, is the sum of 
the electron and phonon momenta, 



K 



q n q 



and it commutes with the Hamiltonian. In the present 
treatment, only the low energy polaron states (the low 
lying states of the system for which the electron and lat- 
tice part of the wave function are spatially bound), are 
explicitly calculated. In this case, the total momentum 
K of the system is also the polaron momentum. 
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III. METHODS 



The eigenstate computations reported here are based 
on the variational approach. Still, from the physical and 
mathematical point of view, there are significant differ- 
ences among them. Each method is therefore described 
separately in the present section. It is important to notice 
that two qualitatively different kinds of polaron functions 
are employed, the localized ones and the translationally 
invariant ones. The primary objective of the methods 
with localized functions are self-trapped polaron states, 
i.e., the strong coupling regime. On the other hand, 
translationally invariant states are devised in the first 
place with the weak coupling and crossover regime in 
mind. The methods with a small number of variational 
parameters are meant to help in understanding the ba- 
sic properties of the polarons in different regimes of the 
Holstein model. Again, the methods with a very large 
number of variational parameters are necessary to obtain 
accurate results for the polaron states. 



A. L method: A simple localized polaron function 

Let us start with a simple localized polaron wave func- 
tion formed as a product of the electron and the lattice 
part, centered at the lattice site j, 

= £ Vn c] + J ([] S 3+m (U) |0> • (4) 

n m 

Here, r\ n is the normalized electron function at site j + n, 
En VnVn = L while Sj+ m (£, m ) denotes the coherent state 
operator acting on site j + m, with a complex amplitude, 

Sj+rnitm) = exp (£ m &j +m - CA+m) • 

It is easy to see that operator Sj+m(£ m ) shifts the space 
and momentum coordinates of lattice vibration at a site 
j + m by 23?(£ m )a:o and 2S(£ m )p , respectively, 



o It \ _ Ji2Q(U)poS:j+ m /K\ [-i2fSt(U)x p j+m /K\ 



x n = (a/ k) g n , 



(6) 



x e 



The variational energy of the state (|J), E v , is indepen- 
dent of j , and is given by 



so that only the equation for r\ n has to be solved. The 
well known approximate solution to this problem is the 
large Holstein polarona valid in the long wave limit, 



r\ n = 



cosh 1 {n/d po i) , dp i = 



2t Hlu 2t 



r-i (7) 

The numerical scheme suggested in Ref. |ll|, and denoted 
here by L (L for localized), is not restricted to long waves, 
and is used here in order to obtain the exact minimum 
of Eq. (H). The energy, henceforth referred to as El, 
depends only on two relevant Hamiltonian parameters, 
t and g 2 /ncu = e p , and therefore El and x 1 ^ are both 
independent of M. For this reason Eq. (||) is sometimes 
referred to as the adiabatic locking of electron and lattice 
coordinates. 

It is noteworthy that the lattice part of the function 
\<fi)f is a simple product of coherent states with real am- 
plitudes. Consequently, the state of the lattice at some 
site is defined by the ground state of the displaced har- 
monic oscillator. The mean lattice deformation cor- 
responds to the equilibrium position of that oscillator, 
while the lattice zero point motion is approximated in 
Eq. (0) by that of the free lattice. 



B. T and CT methods: Translational polaron 
functions and the first excited state 



Next, we shall study a translationally invariant solu- 
tion composed of a linear superposition of the localized 
states (0), 



I* 



K 



(8) 



| 'I' if) describes the polaron state with the momen- 
tum K. A similar type of function was first proposed 
by Toyozawa.tHI In the present work, £ m = g/nuj is 
used so that the mean total deformation of function (||) 
satisfies Eq. (||), 



Xtot = x 



E V = _ *$^ r /n( J 7n-l + J 7n+l) + ^ $Z ' ^™ ^ 
n n 

- gJ2\ r ln\ 2 (tn + Cn)- (5) 

n 

The minimization of the energy with respect to £ n es- 
tablishes a simple relationship between the lattice mean 
deformation and electron density g n = \n n 



l ln 1 2 1 



The expectation value of the polaron energy, E^j, may 
be written in terms of \(fj), 



E„= Ez 



„iKAa 



iKAa 



E 



A ' 



?iKAa 



Sa 



(9) 



Ea and S*a are given in the Appendix. A simple method 
for calculating the minimum of the energy E^, has not 
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yet been proposed, but accurate results have been ob- 
tained in Ref. [l3] by using the Toyozawa method, which 
includes a very large number of variational parameters. 
Some additional approximations may be found in Refs. ^ 
and |l^ , p^|Jl^ , |l7| , p^| . The approximation used here simpli- 
fies the general expression in Eq. (||) by introducing the 
exponential form for functions rj n and £ m ,EJ3 

r?„ = GG 1 " 1 , f m = AB^e lKma , < G, B < 1 . 

_ (10) 
Equation ( |10| ) defines a polaron function l^xiGjB)) 
which is completely determined by two parameters, G 
and B. 

In what follows, two different approaches are pre- 
sented. The first, denoted by the index T (T for transla- 
tional) , treats G and B as the variational parameters for 
which the energy minimum Et has to be found, and its 
corresponding polaron function is l^ 7 ^). The T method 
gives good results in the weak and the strong coupling 
regime. Namely, in both of these limits, the function 
becomes similar to the polaron function obtained 
by the appropriate perturbative calculations. 

In the second approach, presented here for the first 
time, the variational method is used to define a gener- 
alized eigenvalue problem as follows. The polaron wave 
function, denoted by the index CT (CT for combination 
of translational functions), is rewritten as a linear com- 
bination of functions |4 , k(G„, B n )), 

\® C k) =$> n |\MG„,5„)) • (11) 

n=l 

It is understood here that the functions \^(G n , B n )) form 
a set of p generally nonorthogonal functions, defined by 
p different pairs of parameters (G n ,B n ). Again, the co- 
efficients a n should be determined from the requirement 
that the expectation value of the energy, 

{* c K T \m c K T ) 
CT (m* c K T ) ' 

is minimal, 

dE CT /dal = 0, l<n<p, 

or, 

Y J {^K{G n ,B n )\H\^ K {G n ,B n ,)) a n , 

= E CT J2(*K(G n ,B n )\* K (G n ,,B n ,)) a n , . 

n' 

The solution of this generalized eigenvalue problem is a 
set of p orthogonal polaron functions, |3>§>„), with cor- 
responding energies Eqt . The ground state energy is 



Eqj,. One may always include the function l 1 ?^) in the 
sum ([n]) in order to ensure that the energy EjSj, is the 
same or better than Et, the energy computed by the T 
method. Moreover, by paying some further attention to 
the starting set of functions \^>k (G n , B n )) in Eq. ( |Tl| ) at 
the outset, one is able to investigate the first excited state 
\$Km=i) °1 l ne system, when this state corresponds to 
an excited polaron. The best results for the CT method 
are obtained when the number p of \^ K{G n , B n )) func- 
tions in Eq. (|§|) changes with the Hamiltonian parame- 
ters. The special case where the CT method is used with 
constant p = 2 is denoted by the index CT-. 



C. eT and eL methods: Exact translational and 
exact localized polaron functions 

Finally, this paper presents the results fjif.tMJO numer- 
ical exact diagonalization methods .oEEdEj'El In order 
to compute the low energy polaron states, one approxi- 
mates the infinite dimensional Hamiltonian matrix with 
a finite one, and proceeds with the exact diagonalization 
of this matrix. The lowest eigenvalue and eigenvector in 
such a reduced Hilbert space correspond to the polaron 
energy and wave function, respectively. For a large sparse 
matrix, the energy and the wave function can be calcu- 
lated very accurately, by using an appropriate numerical 
scheme, in the present case the Lanczos algorithm. 

The two exact diagonalization methods used here differ 
in the choice of the basis of the Hilbert space. In the first 
method,E3 denoted by eT (e stands for exact diagonal- 
ization and T for translational), the general orthonormal 
state is given by 

|n ,n_i,m, -,n m , ...) e ^ = 
-=^e iJria ct|no ) n_i,ni,...,n m ,...) j , (12) 

3 

which describes an eigenstate of the system with momen- 
tum K. n m is the number of phonons at the mth lattice 
site away from the electron. For example, at the site j, 
and at the nearest neighbor lattice sites left and right 
from it, there are no, and n\ phonons, respectively. 
The Hamiltonian (|l|) does not mix states ( |l2|) with differ- 
ent momenta. Therefore, the polaron function obtained 
by the eT method has the same K momentum as the ba- 
sis states. The current implementation of the eT method 
is highly accurate, and from a practical point of view may 
be treated as exact.E3 For this reason, the eT results can 
be used to determine the numerical errors present in the 
other methods. 

The minimal number of states of the reduced basis nec- 
essary to obtain accurate results depends on the Hamil- 
tonian parameters. This number for the eT method in- 
creases very rapidly for huj <C g, t, which prevents its use 
for both large g and t. 
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In the second exact diagonalization method, denoted 
by eL (e stands for exact diagonalization and L for local- 
ized) , the general orthonormal state of the chosen basis 
is more complicated than for the eT method, 

\i, n ,n_i,ni, ...,n m , £ m )^ L = 

4+t [n^'+m^m)] I n ,n-i,nx,... ,n m , .(13) 

m 

Here the i and m indices are given with respect to the 
center of polaron, which is placed at the site j. Thus cj +i 
creates an electron at the ith site from the polaron center 
at j, n m is the number of extra phonons at the mth site 
away of the polaron center, when the lattice is already 
distorted by the coherent state operators Sj +m (^ m ), i.e., 
Xj+m — 2xo£ m - The eL method calculates localized po- 
laron wave functions. Namely, Eq. ( |l3| ) describes a local- 
ized state, with the polaron center at site j kept constant. 
The eL method is therefore accurate only in the strong 
coupling regime, in which the effects of polaron dereal- 
ization are negligible (self-trapped polarons). 

For a given set of Hamiltonian parameters, £ m in Eq. 
( |l3| ) are determined by the use of the L method, i.e., 
by minimizing the energy (|5|), £ m = £4- If omv those 
states (|l|) with all n m — are used in calculations, the 
eL method gives the same polaron wave function as the 
L method. The additional states dl3|), with n m phonon 
excitations, are necessary to obtain the actual equilib- 
rium positions of the lattice in the exact localized polaron 
state and the zero point motion of the renormalized lat- 
tice vibrations. It is worth noting that the electron and 
phonon parts of the eL polaron wave function cannot be 
completely separated, as in the case of the L function. 

In the case of eL method, the mean lattice defor- 
mation of the localized polaron is approximately taken 
care of by the product of the coherent states opera- 
tors, Sj+m(tm), which keeps the necessary number 
of states (|13| ), in the eL method, relatively small. In 
order to reduce the basis of the Hilbert space, the max- 
imal allowed distance of the electron and phonons from 
the polaron center has been limited here by the choice 
| * |,| m |< D max . The distance D max has been 
determined from the condition that £ m /£o < 10~ 4 if 
m > D max . The maximal total number of phonons has 
been kept limited, retaining the states with ^ TO n m < 4. 
As the sum J2 m n ™ does not include the phonons asso- 
ciated with the coherent state operators S , J+m (^ m ), the 
small value of n m is not a restriction on the overall 
amplitude of the lattice deformation. The accuracy of 
the results obtained by the eL method, supplemented by 
the two abovementioned criteria, depends of course on 
the values of parameters. 

For the purpose of clarity, it seems appropriate at this 
point to review briefly the notation L, T, CT, eT, and 
eL of all five presented methods. All methods based on 
the localized polaron function are denoted by the letter 
L (L and eL methods), whereas all methods based on 



the translational function are denoted by the letter T 
(T, CT, and eT methods). The letter e denotes an exact 
diagonalization method (eT and eL methods), while a 
single letter notation (L and T) suggests the simplest 
form of the method. 



IV. RESULTS 

As the variational methods of the preceding section 
themselves, the results of the corresponding calculations 
may be best understood in terms of the weak and strong 
coupling limits and the-crossover between them. It has 
been shown previouslycJ on the basis of the Global-Local 
method for 0.1 fiuj < t < 10 Tilu, that the empirical rela- 
tion 



9st — fiuo + Vthco (14) 

describes well the values of parameters for which the vari- 
ation of the effective polaron mass with g is the fastest. 
It will be also argued here that gsr of Eq. ( |l4| ) describes 
accurately the crossover from the weak to the strong cou- 
pling limit with respect to the nature of the K = 
ground state. The latter changes continuously from the 
light delocalized state in the weak coupling limit to the 
heavy self-trapped state in the strong coupling limit, with 
the anticrossing of the two states at g ~ gsT- The physi- 
cal content of Eq. (|lj) is best understood by considering 
the limits of small and large t with respect to Ticu, when, 
respectively, gsr ~ and Ep T = g 2 ST /Titu « t. Both 
these conditions were qualitatively explained in Ref. |2^. 
In the present section we first discuss the strong and the 
weak coupling limit, and then devote most of our atten- 
tion to the crossover between the two. 



A. Strong coupling limit 

The nature of the self-trapped polaron may be dis- 
cussed by analyzing the properties of the two polaron 
wave functions given by Eqs. (Q) and (||). When the 
localized polaron functions \ipj) at different lattice sites 
are not orthogonal, 

Sa = (Vj|Vj+A) <5o,A , 

the local properties and derealization effects of 1^^) in 
Eq. (^) are a complex mixture. However, in both the 
weak and strong coupling regimes, the translational po- 
laron function can be approximately written in terms of 
orthogonal localized polaron functions. In the weak cou- 
pling regime, the orthogonality follows from the electron 
part of the wave function, 

n 
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Figure 1: The solid curve is one fourth of the polaron band- 
width, Eq. @. The long-dashed curve corresponds to the 
error of the polaron energy calculated by the eL method, 
EeL^EeT- The short-dashed curve is the electron hopping en- 
ergy reduced by the Debye- Waller factor, Eq. (^). All three 
curves are given as functions of the electron-phonon coupling 
g = gsT + fru = yi hco + 2 hui. The energies are given in units 
of hu (i.e., hco = 1). 



while in the strong coupling regime, it follows from the 
lattice part, 

5 A ~ Y A = cxp [-1 - £ m+A ) 2 ] « 6 0A ■ (15) 

m 

In Eq. (|l5|), Ya is the Debye- Waller factor. The condi- 
tion ( filf ) corresponds to the regime of self-trapped po- 
larons. The negligible contribution of the lattice part 
to the overlap of any two localized polaron functions at 
different lattice sites results then in a negligible polaron 
hopping energy t po i- 

One may notice that in the limit Ya^o - ► the trans- 
lational form of the polaron function (||) has no conse- 
quences on polaron energy, and the minimal values of 
the variational energies (||) and (^|) coincide, i.e., Ej\, 
Sa ~ ^o,A m Eq. (||). The hopping of the self-trapped 
polaron occurs only at the time scale which is much larger 
than the scale relevant for the local interplay between the 
electron and the lattice deformation. Therefore an accu- 
rate description of the local polaron properties may be 
obtained even if the hopping of the polaron is completely 
omitted by using only localized polaron functions. The 
translational invariance of the polaron may be, however, 
always restored in the same way as the function (^) is 
used to obtain (||). 

The eL method provides very accurate results for the 
ground state energy of the self-trapped polaron. Ex- 
cept that it neglects the polaron hopping, the local po- 
laron function is calculated exactly. The error of the eL 
method may be estimated by using the eT method, that 



is, by subtracting the exact energy of zero momentum po- 
laron state E e T from E e L . For the electron-phonon cou- 
pling which is greater than the critical electron-phonon 
coupling in Eq. (|l4| ) just by hu, i.e., g — gsr + hu), 
the maximal error of the eL method is E e L — E e x < 
3 x 10~ 4 huj. This is shown in Fig. [l], in which two quali- 
tative estimates of t po i in the strong coupling regime are 
plotted as well. The first estimate gives t po i as one fourth 
of the polaron bandwidth computed exactly by the eT 
method, 



pol 



-W 



eT 



■[E eT (K = n/a)-E eT (K = 0)], (16) 



while the second estimate, based on the L method, mul- 
tiplies the electron hopping energy with the small Debye- 
Waller factor, 



(17) 



tpoi ~tY 1 L =t CXp [-1 " ^n+lf 



The estimation of the Debye- Waller factor in Eq. (|17]) is 
based on the evaluation of the lattice part of the over- 
lap of the two neighboring \<p^) functions (the L method 
gives good results for the mean lattice deformation of 
the self-trapped polarons). All three curves in Fig. [l]are 
similar functions of g, which is not surprising since all 
plotted quantities are related to the polaron hopping en- 
ergy t po i in the strong coupling regime. Moreover, it may 
be seen from Fig. [j] that the error of the eL method is 
almost equal to one half of the polaron bandwidth, 

E eL - E eT w W eT /2 . 

One of the advantages of the eL method is that it per- 
mits separate calculations of the electron and the lattice 
properties, in spite of the fact that the electron and the 
phonon part of the eL wave function are not separable. 
For instance, for the polaron centered at the origin, the 
associated mean lattice deformation x n is given simply 
by the expectation value of x n of Eq. (|^). Besides x n , in 
the present paper the mean real space uncertainty of the 
on-site lattice vibration Ax n and the product Ax n Ap n , 
where Ap n is the mean momentum uncertainty of the 
on-site lattice vibration, arc calculated. 

Figure || shows the data for two sets of parameter, 
corresponding to a small self-trapped polaron and a self- 
trapped polaron extended over few lattice sites, respec- 
tively. For the second set of parameter (large g and t), 
the eT method is not accurate enough, and the question 
arises of whether the eL polaron state is really a good 
approximation of the system ground state. It is difficult 
to prove that a polaron state, with non-negligible hop- 
ping energy t po i and an energy close or lower than eL 
energy, does not exist. This question is currently under 
investigation. 

Differences between the results of the L and the eL 
methods are also analyzed. It is found that they are more 
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Figure 2: Difference between the mean lattice deformation, mean uncertainty of the on-site lattice vibration and corresponding 
product of uncertainties for the L and eL methods. The inset shows the mean lattice deformation of the eL method. Hamiltonian 
parameters are t = 10 Tiui, g = 4.5 hco, and t — 250 hcu, g = 16.5 Tiuj, for the first and second plot, respectively. 



pronounced for small g. We may see from the results in 
Fig. |2|, that the mean lattice deformation differs between 
these two methods. In the case of the eL method, the 
mean lattice deformation is more extended and the width 
of the polaron is slightly larger. Additionally, the elec- 
tron density of the eL method remains approximately 
proportional to the mean lattice deformation (as in Eq. 
(||), valid for the L method), since 

K L - ^Q* L Xo]/< L < 1% . 

As has already been mentioned, the lattice part of the L 
function describes a set of displaced, but unrenormalized, 
harmonic oscillators, so Ax„ — xq, Ap^ = po, where Xq 
and po are defined in Eq. <^). One may notice from Fig. || 
that real space uncertainties of the lattice vibrations on 
the sites occupied by the polaron are larger in the eL case 
than in the L case, Ax^ > x , but uncertainties of on- 
site momentum lattice vibrations are smaller, Apf^ < pq. 

It may be concluded, from the first plot in Fig. [| that 
the electron affects mostly the central site of the small 
polaron. The product of uncertainties for this site stays 
close to the free lattice value, Ax^l Ap^i w xoPo- 
Therefore the phonon mode at the central site of the 
small polaron may be treated, in a good approximation, 
as harmonic. In particular, the renormalized frequency 
of this mode, wf t £ ' can be roughly estimated from the 
relation 

A<£ a • 

Since Axf^± Q > Xq, it follows that w^£ < w - Therefore in 
the strong coupling regime the first excited state should 
correspond to the excitation of the renormalized phonon 



mode, rather than to the excitation of the phonon of en- 
ergy hu>, which is uncorrelated to the polaron. One may 
also notice that the energy of the mean lattice deforma- 
tion is larger for the eL method than the L method (for 
the L method this energy is minimal). This is compen- 
sated for, however, with the lower energy associated with 
the zero point motion of the Cj^Lq phonon mode, which 
makes the total polaron energy of the eL method lower. 

For the more extended self-trapped polarons, the 
renormalized normal phonon modes are expected to be 
spread over a number of lattice sites. Consequently, a 
number of different phonon modes contribute to the lat- 
tice displacement at the lattice sites occupied by the po- 
laron. Thus the analysis of the on-site vibrations can- 
not give direct information on the renormalized lattice 
modes. Nevertheless, from the second plot in Fig. || one 
may notice that the product of uncertainties Ax^Apf^ 
shows a minor deviation from that of the harmonic vibra- 
tion. This suggests that the renormalized lattice modes 
of the extended self-trapped polaron are harmonic as well 
in a good approximation. 



B. Weak coupling regime 

In the weak coupling regime, the T method gives re- 
sults close to the eT results. Since the form of the func- 
tion T is quite simple, it will be used in this section as a 
basis for further discussion. In the weak coupling regime, 
the minimum of the variational energy Et corresponds 
to small values of the variational parameters G and A. 
The standard perturbative ground state of the system 
with momentum K in terms of the T method polaron 
function may be written as follows: 
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(l+Aj2B lm{ e iKma b} + J |0) , if < X c ,(18) 



x(l + /l^Bl™li,t +iJ | 0)> JT> Jf c . (19) 



For K < K c the wave function ( |l8| ) has two parts. The 
main part corresponds to the free electron of momentum 
K, and the smaller part, proportional to A, corresponds 
to the electron dressed by one spatially correlated vir- 
tual phonon. At the threshold K c , the energy of such 
a polaron state intersects with the energy of the system 
consisting of the zero momentum polaron and one extra 
phonon with momentum K c , Eq. (p9|). So, for K > K c 
the ground state is achieved with one real phonon in the 
system which carries the system momentum and which is 
spatially uncorrelated with the polaronlij'E] For K < K c 
this state becomes the first excited state of the system. 
The difference between the energies of the ground state 
and the first excited state is the largest for K = 0, and 
is equal to fko. 

The validity of the perturbative treatment requires 
that the weight of the second term in Eqs. ([l8]) and ( |l9| ) is 
small, i.e., the mean number of phonons associated with 
the lattice deformation has to satisfy 



N' 



It = A 2 (l + B 2 )/(1-B 2 ) 
= , 9 2 (l-j?)(l + 5 2 ) 

(M 2 (i + s) 3 



< 1 . 



(20) 



Here, A has been eliminated by using Eq. (g). There are 
two ways to satisfy condition (Ej). Either the electron- 
phonon coupling is small, g <C tkJ, or the lattice de- 
formation is spread to a large number of lattice sites, 
1 — B <C 1. In the latter case, the total mean polaron 
deformation does not have to be small, x to t/xo — 2g/hui, 
since g can be larger than Huj. 

The translationally invariant form of the wave function 
for the T method, given by Eq. (Q), provides an energy 
gain due to the polaron derealization. At the same time, 
the spatial correlation between the electron and the lat- 
tice deformation has a finite length. For instance, in Eq. 
(||), this length is of the same order for \^k) and \<pj). 
The perturbative calculation for B in Eq. (Jlq) gives 



B = cos (Ka) + hu/2t - y/[cos(Ka) + hw/2t} 2 - 1 . 

B measures the electron-lattice correlation length. B is 
independent of g, which makes the correlation length fi- 
nite, even in the limit g — > for which the lattice defor- 
mation vanishes, A ~ g/tku. 



On the other hand, the electron-lattice deformation 
correlation length and the polaron derealization range 
are of the same order for the localized functions. This 
may be easily seen from Eq. (||). In the limit g — ► they 
both become infinite. Thus in the weak coupling regime 
one obtains a localized polaron state \<pj) of very large 
width, but associated with a tiny lattice deformation. 
This is specific for one dimensional systems in which an 
attractive symmetric potential has always a bound elec- 
tron state. Therefore the corresponding polaron energy is 
less than the free electron energy — 2t. In higher dimen- 
sions, an arbitrary attractive symmetric potential does 
not have a bound electron state for sufficiently small g 
(or the electron binding energy is too small to balance the 
lattice deformation energy), and the total polaron energy 
is larger than —2dt, where d is the dimension of the sys- 
tem. This explains why, in the weak coupl ing r egime, 
the adiabatic (localized) polaron functionsO'EZlca fail to 
have an energy lower than the free electron energy, for 
the dimension of the system greater than one. A detailed 
parallel perturbative analysis of polarons in one, two, and 
three dimensions is given in Ref. B9. 



C. Crossover regime 

In the crossover regime the polaron hopping energy 
tp i is not negligible, and the translationally invariant 
polaron functions should be used in order to obtain the 
full physical picture of the polaron. In order to calculate 
the numerical errors of different methods accurately, the 
present discussion of the crossover regime is restricted to 
the values of t smaller than 25 hu>, for which results from 
the eT method are available. 

In order to examine the crossover regime it is instruc- 
tive to calculate the energy difference between the ground 



and first excited state, Eqj, — Eqj,. Let g c denote the 
value of g for which this difference is minimal, 



(0) 



d{E^ T -E^ T )/dg\ 







Our results indicate that g c is very close to the value 
of gsr given by Eq. For example, for t = 20 hi*), 

g c = 5.55 htv, while gsr = 5.47 Tiuj. For smaller t, g c and 
g ST coincide even better. The analysis of the effective. 
mass,E-3 variational energy of the polaron ground state,E2l 
polaron size,Ed as well as the behavior of the first excita- 
tion energy, suggest that a dramatic change in the nature 
of the polaron ground state and the first excited state are 
intimately related in the crossover regime. 

This can be well understood by considering the prop- 
erties of the wave function in the T method. Near g Cl 
this polaron wave function has two separate energy min- 
ima in the G-B parameter space defined by Eq. (px|), 
which become degenerated for g = g c . Let the symbol 
< denotes the lower minimum at g < g Cl and the sym- 
bol > the lower minimum at g > g c . and |^ ,> ) 
are the corresponding polaron wave functions. Even if 
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3.55 g 3.75 



Figure 3: The anticrossing of ) and |$ r> ) energies together 

with the ground state |$ _ ) and first excited state |<E>i _ ) 
energies are shown. The results of the CT method are also 
plotted for comparison. The ground state energy obtained by 
the eT method is subtracted from all energy curves, t = 5 Tloj, 
K = 0, and Tilu — 1, while gsr = 3.24 Ttw. 



they are not mutually orthogonal, they are still physi- 
cally quite different. The numerical data show that the 
transnational invariance of |^' < ) contributes strongly to 
the polaron energy. On the other hand, the translational 
invariance of | 1 3/ > ) has almost negligible energy contribu- 
tions, i.e., \^f > ) describes an almost self-trapped polaron. 
It is worth noting that the degenerate nature of the vari- 
ational eaergy which has been reported for the Toyozawa 
methocO is of the same kind as the one of the T method 
discussed here. Namely, both T and Toyozawa method 
are based on the same polaron function (g). 

| X I ,< ) and \' i & > } can be combined to form new polaron 
functions, 



I* 



CT- 



(21) 



It should be mentioned that a similar combination of two 
states has been already used in Ref. |3p in order to cal- 
culate the polaron ground state. |$ CT_ ) corresponds to 
Eq. ( pd| ) with p = 2, which means that the CT- method 
is implied. From this treatment the improved ground 
state \&q T ~) and the approximate first excited polaron 
state |$f T_ ) are obtained. It may be seen from Fig. || 
that CT- method describes an anticrossing of \^f > ) and 
|\E ,< ) states, which yields two orthogonal states, \^q T ~) 
and \&f T ~). In order to have a better illustration of that 
anticrossing in Fig. ||, the exact ground state energy E e x 
is subtracted from all plotted energy curves. For g < g c 
the light state \^ < ) is lower in energy, and participates in 
the ground state more than the heavy state | X I' > ). How- 
ever, as g increases this balance changes continuously in 



favor of | \E' > ) . The opposite trend is observed for the first 
excited state, which is heavier than the ground state for 
smaller g and lighter for larger g. 

The CT method gives better results for the ground and 
first excited state when a large number of functions (large 
p) in Eq. ( |i"i"| ) is used. From Fig. |3[ one may estimate that 

E^l - ~ hu/2 at g = g c ™™ 



E$l, unlike E { ^ Tl 



satisfies E, 



(i) 

CT 



Moreover, the energy 
< + Huj for all 



g > g c , which is an important improvement over the 
p = 2 result. 



V. SUMMARY OF RESULTS 

Figure ^ shows the results for the ground state energy 
and the energy of the first excited state with total system 
momentum K = 0, as functions of g, obtained by several 
different methods. Two plots correspond to t = 5 hto 
and t = 10 hu, respectively. The ground state energy 
obtained by the eT method is subtracted from all the 
other results. <?i and g2 are used to mark three differ- 
ent polaron regimes with respect to the strength of the 
electron-phonon coupling. For g < gi the mean number 
of phonons of the lattice deformation is less than one. 
Thus for g < g\ we recognize the weak coupling regime. 
For g > g2 — gsr + Tilj the eL polaron energy has a 
negligible error (see Fig. |l|) , which means that for g > g^ 



polarons are self-trapped, and we recognize the strong 
coupling regime. The crossover regime is found in the 
interval gi < g < g 2 , with gi < g S r ~ g c < 92- 

It may be noted from Fig. ^ that in the weak coupling 
regime, the energy obtained by the L method (£"l) is 
close to the free electron energy — 2t (for t > 20 hu, the 
absolute error of the L method becomes greater than hu 
in some intervals of g). In the strong coupling regime, 
for large <?, El approaches the exact polaron energy. 

The error of the T method, as may be seen from Fig. |^, 
is the largest in the crossover regime in which the results 
can be improved by a better choice of r\ n and £ m in Eq. 
(|io|). In the strong coupling regime, the translationally 
invariant form of the function T has no effect on the 
polaron energy, i.e., both the T and L methods give very 
similar results. 

The energy of the first excited polaron state, 
intersects the energy of the ground state plus one phonon, 

E e T + hto, for g > gi. After the minimum of Eq^, — Eq^ 

is reached in the crossover regime at g — g c w gsT> ^ct 
approaches E e T + hu asymptotically for g > g 2 . As has 
already been pointed out in Sec. IV, the excitation of the 



renormalized phonon mode explains the nature of the 
first excited state in the strong coupling limit. In Ref. |28| 
perturbation theory was used to calculate the frequency 
of this local phonon mode, Hi' , to the lowest order in 
t hco/g 2 , 



Co 1 = uj s/l - (t hco/g 2 ) 2 



(22) 
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Figure 4: The ground state energy of the polaron for various methods and the first excited state energy of the CT method are 
plotted for t = 5 hui and t = 10 hui, respectively. K = and hu — 1 for both plots. The eT ground state energy is subtracted 
from all results. Only the lowest hui energy interval of the spectrum, relevant for the ground and first excited energy, is shown. 
The E e T + hui line denotes the first excited state energy, when it consists of the polaron ground state and one extra phonon: 
see Eq. ©. 



The obtained perturbative correction to the phonon fre- 
quency is adiabatic, i.e., the square root in Eq. (|2^) is 
independent of mass M . It is worth noting that for large 

g, E^,)p calculated here shows approximately the same 
behavior as the energy of the ground state with one ad- 
ditional phonon of frequency ( |22] ) , 

« E eT + hCo' . 

The lattice part of the CT function is spatially sym- 
metric with respect to the electron. Therefore in the 
strong coupling regime the first excited state, which has 
been identified here as a local renormalized phonon mode 
of the self-trapped polaron, should be basically a sym- 
metric oscillation of the lattice deformation around the 
central polaron site. In the vicinity of gsT, on the other 
hand, the nature of the first excited state is currently 
explained by the anticrossing of the self-trapped and the 
delocalized polaron state. The question of how the ex- 
cited self-trapped polaron state in the strong coupling 
regime and the first excited state near gsr may be linked 
together is a matter of further considerations. One may 
speculate that the anticrossing of the excited self-trapped 
polaron state and the delocalized state will give the an- 
swer. 

VI. CONCLUSIONS 

The present paper discusses the ground and first ex- 
cited states of the polaron for three different regimes of 
the electron-phonon coupling parameter g. The results 
can be briefly summarized as follows. In the strong cou- 
pling regime the polaron hopping energy to neighboring 



sites is negligible, and the self-trapped polaron states are 
obtained. The results of the eL method suggest that the 
adiabatic picture of the localized polaron state is valid, in 
which some of the local lattice vibrations are renormal- 
ized by the presence of the electron. The numerical data 
show no significant deviation from the adiabatic locking 
relation (^) of the electron site density and the mean lat- 
tice deformation. In the small polaron case, the predomi- 
nant effect of the electron is the lowering of the frequency 
of the vibration at the central polaron site. The excita- 
tion of the renormalized phonon mode corresponds to the 
first excited state of the small self-trapped polaron. 

The nature of the polaron ground state in the crossover 
regime has been discussed in a number of papers, and its 
rapid change with g has been well established numeri- 
cally. The difference between the energies of the first 
excited state and the ground state, as a function of g, 
has a minimum for g — g c . It is shown, by using the 
CT- method, that the anticrossing of the self-trapped 
and the delocalized polaron state can link the behavior 
of the ground and first excited polaron state. Accord- 
ing to the CT method, for g > g c the effective mass of 
the ground state is larger than the effective mass of the 
first excited state, while for g < g c the opposite is true. 
In addition, it is found that g c , which characterizes the 
first excited state, and gsT, obtained from ground state 
analysis (Eq. (|l4|)) almost coincide for t < 25 hu. 

Upon further reduction of g the total mean number 
of phonons bound by the polaron becomes smaller than 
one, and the weak coupling regime is reached. The nearly 
free electron is dressed by a cloud of virtual phonons, and 
its mass is slightly renormalized. The first excited state 
of the system, with momentum K « 0, can be viewed 
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as the ground state of the polaron plus one additional 
uncorrelated phonon, rather than as an excited polaron. 

Finally, it is worth noting that there is a simple sum 
rule for the mean total lattice deformation, Eq. (Q) , which 
is valid for any number of electrons and is independent 
of the system dimension. This important sum rule may 
be extended to some other models in which the lattice 
deformation is linearly coupled to the local electron den- 
sity. 
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Sa = Ya ^ vlVn+A , 
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Ea = —t Ya ^ VniVn+A+l + Vn+A-l) 

71 

+ huj S A C£ 

m+A 

i n 

- gY A ^VnVn+AiC + Cn+A) ■ 



Appendix A 

In order to calculate 5a and Ea of Eq. 
ing expressions may be used: 



the follow- 
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